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Abstract 

We consider the Friedel sum rule in the context of the scattering theory for the Schrodinger 
operator — + V(x) on graphs made of one-dimensional wires connected to external leads. 
We generalize the Smith formula for graphs. We give several examples of graphs where the 
state counting method given by the Friedel sum rule is not working. The reason for the 
failure of the Friedel sum rule to count the states is the existence of states localized in the 
graph and not coupled to the leads, which occurs if the spectrum is degenerate and the 
number of leads too small. 

PACS : 03.65.Nk, 72.10.Bg, 73.23.-b 

1 Introduction 

This article follows jl| in which we have considered the scattering problem for the Schrodinger 
operator on graphs. The graphs we are interested in are networks made of one-dimensional wires 
identified with finite intervals of M and being connected at vertices. The study of such systems 
has been shown to be relevant in many contexts (for references see |2|, |3], [|, & [f| [7|, [lj). For 
example graphs have been often considered as simple modelizations for the mesoscopic networks 
realized experimentally. In this context the scattering theory is a fundamental tool involved in 
the study of transport properties and many other questions. Several works have been devoted to 
the study of scattering theory on graphs among which we can quote ]|, ||, ||, [| . In our work we 
examine an important aspect of scattering theory, namely the relation between the scattering 
and spectral properties that is established through the Friedel sum rule (FSR). The essence of 
the FSR is to count the number of states in the scattering region, that is related to the phases 
of the eigenvalues of the scattering matrix. The purpose of this article is to show that one must 
be careful when applying the Friedel sum rule to graphs since this formula does not hold for any 
graph, one of the reasons of the breakdown being the occurence of degeneracies in the spectrum 
of the graph (this is a necessary but not a sufficient condition as we will see). 
Here we consider the Schrodinger operator 

H = -D 2 X + V(x) (1) 

where D x = d^ — \A{x) is the covariant derivative ; the x coordinate lives on the graph Q. We 
briefly recall the notations adopted in previous works @, [l(], |l). The graph Q is made of B 
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one-dimensional wires connected at V vertices. We will denote the vertices with greek letters 
(a, /?, n,. . . ). The V x F-adjacency matrix a a p (or connectivity matrix) is defined as: a a p = 1 
if the vertices a and (3 are connected by a bond and a a p = otherwise. The coordination of 
the vertex a (number of bonds issuing from the vertex) is therefore m a = Y^f3 a a/3- The bond 
between the vertices a and /3 will be designated with parenthesis: (a/3). We also introduce the 
notion of arc which is an oriented bond. Each bond (a/3) is associated with two arcs: a/3 and 
/3a. The arcs are labelled with roman letters (i, j,. . . ) and we denote the reversed arc of i with 
a bar: i. 

The coordinate x a p on the bond (a/3) of length l a p belongs to the interval: x a p 6 [0, l a p] 
where l a p is the length of the bond (a/3) (note that by definition xp a = l a p — x a p). 

The Schrodinger operator acts on scalar functions ip(x) living on Q that are represented by 
a set of B components ip[ a i3){ x ap) satisfying appropriate boundary conditions at the vertices 
§0: 

(i) continuity 

VVft) {xafc = 0) = ip a for i = 1, • • • , m a (2) 

{Pi j i = 1, • • • ,m a } is the set of vertices neighbours of the vertex a ; the wave function at the 
vertex is ip a . 

(ii) A second condition sufficient to ensure current conservation (i.e. unitarity of the scattering 
matrix) 

^2 Ga P Dx <*f3 ^(«/3) ( X «/3 = 0) = Klpa , (3) 
& 

where X a is a real parameter. Due to the presence of the connectivity matrix a a p, the sum runs 
over all neighbouring vertices linked with vertex a. 

Note that the conservation of the current alone leads to more general boundary conditions 
and does not require the continuity of the wave function at the nodes (see |l2| for example) 

The magnetic flux along the bond is denoted by 9 a p = f^dxA(x) = —0p a . 

In a scattering situation the graph is connected to the external by leads plugged on vertices 
(we designate by "graph" the compact part that does not include the leads). We call L the 
number of leads through which some plane wave is injected. The quantity of interest is the 
on-shell scattering matrix S which is a L x L matrix that relates the incoming amplitudes in the 
L channels to the outcoming ones. We call (resp. -B„ xt ) the incoming (resp. outcoming) 
amplitude on the external lead connected at the vertex a (i.e. is the coefficient of a plane 
wave e~ lkx sent from the lead connected to vertex a). By definition: 

B cxt = S A cxt . (4) 

The purpose of [|| was to formulate in general terms the scattering theory for the Schrodinger 
operator, generalizing results known in the absence of potential V(x) (Laplace operator) P, pi 



13]. We have found various expressions of S for arbitrary graphs and related S to matrices 
encoding informations about the topology of the graph, the potential on the bonds and the way 
the graph is connected to leads. We recall here the main results of that will be necessary in 
the following discussion. 

(A) The arc matrices formulation. 

We express here E on the energy shell E = k 2 in terms of matrices that couple arcs. The graph 
is described by 2B internal arcs. L external arcs are associated to the L leads. We introduce 
the matrix R that encodes the information about the potential on the graph and couples the 
2B internal arcs of the graph: 

Rij = n Sij + ti 6- ij:j (5) 
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is the matrix element between arcs % and j. The potential on each bond (i) is characterized by 
reflection and transmission coefficients: rj, ti for the injection of the wave in the direction of arc 
the i and rj, t\ for the injection in the direction of the reversed arc i. R is the bond scattering 
matrix. If the potential vanishes (V(x) = 0) we have = and ti = exp(ikli + i9i). 

Next we introduce the vertex scattering matrix Q that encodes the information on the 
topology of the graph and the way it is connected to leads: 

2 

Qa = — — - — 1 if i = j (i issues from the vertex a) (6) 

m a + iX a /k 

2 

= — — if i ^ j both issuing from the vertex a (7) 

m a + \\ a /k 

= otherwise . (8) 

This expression of the vertex scattering matrix is a consequence of the conditions We 
have also explained in |l| how this matrix is affected by the introduction of tunable couplings 
to the leads. The {2B + L) x {2B + L)-matrix Q couples the 2B internal arcs together but also 
the latter to the L external arcs. If we separate Q into corresponding blocks: 



Q 

then the scattering matrix is: 



Q' mt 


Q T 


Q 


Q ext 



(9) 



S = Q ext + Q(R ] - Q int ) _1 Q T . (10) 

The expression ( |l0| ) generalizes the result known in the absence of potential 
(B) The vertex matrices formulation. 

The previous approach is quite natural since we considered scattering matrices of the different 
parts of the system but it has the disadvantage of dealing with rather big matrices. It is more 
efficient to consider matrices that couple vertices. We define the L x ^-matrix W that encodes 
the information about the way the graph is connected to leads: 

W a /3 = w a 5 a /3 (11) 

with a £ V cx t and (3 £ V, where V = {1,---,V} is the set of vertices and V ex t the set of 
vertices connected to leads (Card(V ex t) = L). The parameter w a S M describes the coupling 
between the graph and the lead at vertex a ; its precise physical meaning is discussed in Q. 
In the arc matrices formulation these parameters are introduced in the matrix Q [Q]. We just 
recall that w a = 1 corresponds to perfect coupling [the case considered above in paragraph (A)], 
whereas w a = corresponds to disconnect the lead. The limit w a = ±oo also corresponds to 
disconnection of the lead, however the current is not allowed to flow through the vertex in this 
case and this way to disconnect the lead is equivalent to impose a Dirichlet boundary at the 
vertex (A a = oo). 

We also introduce the matrix M that contains all the information on the isolated graph 
(potential on the bond and topology): 



M, 




a/3 — "a/3 "t" 7 



(1 ^afi) (1 Tfia) tafj, tfj,a 



2 t a j3 

^al3~? w \ • (12) 

(1 + r a p){l + r Pa ) - t a/3 tp a 
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Then, the scattering matrix reads: 

E = -1 + 2W (M + W T W)~ 1 W T . (13) 

These equations generalize the result known in the absence of the potential [^, |j] . In this latter 
case we recover from (|lj) the well-known matrix: 

i e lda P 

M a /3 = i <5 Q/ 3 2^ °V Cot S kl a^ ~ «a/3 • , , • (I 4 ) 

Some examples of application of these formulae are given in [|J . 

We describe the organization of the paper. Since there has been recently some confusion in 
the literature about the content of the Priedel Sum rule, we think it is useful to spend some time 
by reviewing some aspects around this relation, which will be also necessary for the following. 
In section || we generalize the Smith formula for graphs. Then we provide in section || several 
examples of violation of the Friedel sum rule and explain the origin of this failure. 



2 The Friedel Sum Rule 



To be precise we consider the scattering theory for the Schrodinger equation on a three dimen- 
sional Euclidean manifold and restrict ourselves to the case of a rotational invariant potential 
supposed to be concentrated in a sphere of radius R. A basis of eigenstates is given by the partial 
waves ipi(r)Y l m (9,ip) [where Yi n (9,ip) are the spherical harmonics] whose radial parts involve 
the phase shifts rji(E): ipi(r) = £ sin(Ax — ln/2 + r]i) for r ^ R. The energy of this eigenstate 

is E = k 2 . The Krein- Friedel relation relates the variation of the density of states (DoS) to scat- 
tering properties. We introduce the local density of states (LDoS) p(f;E) = (f\5(E — H)\f). 
We denote po(r; E) the LDoS in the absence of the potential. The relation reads: 



/ 



dr [pif- E) - p (r; E)] = ~ f> + 1)^| . (15) 



Since the integral in the l.h.s. runs over the whole space, the total density of states (DoS) is 
diverging like the volume of integration; however the difference of the l.h.s. is a finite quantity. 
The demonstration of ( |l5| ) in the one-dimensional case is recalled in appendix [B[ 

Using the fact that e 2l7?i are the eigenvalues of the on-shell scattering matrix £ we can writeQ 

J dr[p(r;E)-p (f;E)} = ±±Tr{ln±(E)} = J- A hidet ±{E) , (16) 

where the trace is computed on the energy shell E over channel indices. 

It is convenient to introduce the Friedel phase defined as: 8* (E) = — ilndet S(-E) with the 
additional constraint to be a continuous function of the energy. It is the sum of the cumulative 
phases of the eigenvalues e 1<5a of the scattering matrix S: 8? (E) = ^ a 8 a {E). The Friedel phase 



1 The scattering matrix E entering into (|L6|) is slightly different from the scattering matrix E we have intro- 
duced, both being related through a simple transformation. The phase shifts r/i(E) (phases of the eigenvalues e 2 "" 
of E) encode the effect of the scattering potential compared with the free case: in the absence of the potential, 
the phase shifts r\i vanish. On the other hand the phases 5i(E) of the eigenvalues e" 5 ' of E are measured from the 
edge of the scattering region: ipi(r) = i sin(fc(r — R) — ln/2 + <5;/2) for r ^ R. 

Therefore we have: Si — 2rji + 2kR. The relation between E and E is also explained in detail in the one-dimensional 
case in appendix |b|. 
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counts the number of resonance peaks: if they are sufficiently narrow to be well separated, in 
the neighbourhood of a resonance, the determinant behaves like 



isf(E) • - ■ fE — E n — iA, 



e-^=det£(£) oc ^ 7 , (17) 

up to a constant phase. £7 n is the position of the resonance, A n its width and d n the degeneracy 
of the state. This expression shows that the phase 5^ (E) makes a jump of 2 r Kd n when E crosses 
the resonance. 

The relation (|l^) was derived long ago by Beth and Uhlenbeck |T3 in the context of the 
study of a gas of interacting particles, where it is involved in the second virial coefficient (related 
to the two-body problem). The generalization for a systematic expansion of the grand potential 



was provided in [la]. The demonstration of the Krein-Friedel relation [16, 17, |18L [19[, also called 
Friedel sum rule, is given in standard textbooks for rotational invariant potentials [2^, |2lj|. It 
is also worth mentionning the existence of a vast literature in mathematical physics dealing 
with the scattering theory. Many references can be found in j2^] which devotes its last chapter 
to the study of the Krein spectral shift function (the Friedel phase) and trace formula. The 
matrix -iStg whose trace is computed in ( p~6| ) is the matrix of Wigner time delays (note also 
the existence of a classical formulation of the second virial coefficient in terms of classical time 



delays in pq]). It is worth mentionning that fllq ) is exact which is the beauty of this relation 
(its validity is not restricted to a high energy regime for example). Integrated over the interval 
of energy below the Fermi energy, give the accumulation of charge due to the presence 



of the potential, to use the language of [16]. 

Instead of considering the variation of the DoS of the whole space, it is also possible to study 
the LDoS integrated over the interacting region only. This quantity is also related to scattering 
properties through the Smith formula |24] which defines the time delay. For a rotational invariant 
potential the relation reads: 

2n j\rr*\Mr)\ 2 = + | - -L S m(2kR + 2 Vl - Ztt) = ^| - -L sin (^ - In) . (18) 

Note that this relation was also derived in [16] as an intermediate result for the demonstration 
of (Jl5|) . With a summation over the angular quantum numbers^, we get the LDoS integrated 
over the sphere: 

drp(r;E) = 5>* + l) / drr 2 |^(r)| 2 = - ^(2/ + l) - — sin(^ - /vr) . (19) 
<R l=Q Jo n l=0 \ / 

If the coupling between the scattering region (sphere of radius R) and the external is adjustable, 
this quantity corresponds to the DoS of the scattering region when it is isolated. If we are 
interested in the Weyl contribution of the DoS of the scattering region we can forget the second 
term of ( |l9| ) and consider only the Weyl term of the Friedel phase. 

Due to the central position of the scattering approach in mesoscopic physics, the FSR plays 
an important role in the study of many physical quantities: for example the FSR allows to 
relate the persistent current to scattering properties |2j| and is also involved in electrochemical 
capacitance [26], |27]]. A local formulation was also developed in pq, E3] to relate the local density 
of states to scattering properties (a general discussion of the role of the local FSR is provided 
in p9|). Since graphs are widely used to model mesoscopic networks they have been considered 



2 The choice of normalization for the stationary scattering states ifiE,i,m = 'M'')}'?™^, V 3 ) introduced above 
corresponds to associate to those states a measure dE [it implies for example that J drtp%i m t[>E',l',m' = 
5i lV 5 m , ml 8{E-E')]. 
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to apply concepts involving the FSR, like in K| for the persistent current in a loop connected 
to one lead, or in the recent work |3(]] in which graphs provided examples to illustrate a general 
discussion about some subbtle point related to phases. 

Recently there has been some confusion about the FSR in [31]. Starting from a misinterpre- 
tation of the FSR, these authors claimed that the relation does not hold in the one-dimensional 
case if the potential is made of two 5 peaks, which is not true. We repeat that the general 



demonstration of ( |i6[) in [15] covers the one-dimensional situation. The one-dimensional case is 
reviewed in detail in appendix |B| where we consider as an example the case of one 5 peak (it is 
not difficult to check that the FSR works perfectly well, as it should, for two 5 peaks, a little 
exercice following the same lines). 

The FSR has been proven in arbitrary dimension however it has not been demonstrated for 
graphs which are intermediate objects between one dimensional and higher dimensional systems. 
Then it is important to clarify some points in this context. The FSR (|l5|,16) counts the variation 
of the DoS due to a scattering potential. The Smith relation (^) measures the LDoS integrated 
in the scattering region. Both are based on the idea to count the number of states of the 
scattering region by counting the resonance peaks of the phase shifts derivatives. We will show 
that this procedure is not always applicable for graphs: for example in the case of the complete 
graph that will be studied in detail below, some states of the isolated graph are not manifesting 
by a resonance peak or give rise to a resonance peak that does not carry the correct spectral 
weight (the degeneracy of the level). To study this problem it will be sufficient to consider the 
Weyl part of the Friedel phase to notice some discrepancy with the Weyl part of the DoS of the 
graph. Before following this program and to settle the discussion on more precise grounds, we 
will generalize the Smith formula (18) to the case of graphs, despite it does not always concern 
the DoS as we will see. 



3 Generalization of the Smith relation 



The Smith relation was derived for a one-dimensional system with one scattering channel [24 
(or rotational invariant potentials in 3 dimensions) and involves the Wigner time delay p2| . We 
generalize this relation to the case of a perfectly connected graph (w a = 1). The starting point 
is to introduce 

which satisfies the following relation: 

±n(x) = \m\ 2 (2i) 

dx 

for any solution ip of the Schrodinger equation (— D^ + V(x))ip(x) = Eip(x) ; we recall that 
D x = da — iA(x) is the covariant derivative. Applied to a graph, the relation ([H]) should be 
written for the B components of the wave function on the bonds (and also on the L leads). 

We first derive two properties involving ^l a f3{x a p), the quantity (|20|) related to the component 
of the wave function ijj( a p)(x a /3) (we will denote Oi eac j ^(x) the one associated with the component 
'i/'lead fi (x) on the lead attached to the vertex fi). Obviously, we have £l a j3{x a p) = —£lp a {xp a ). 
Due to the conservation of the current, ensured by (||), the sum of all r^/j's associated to the 
arcs issuing from the vertex /j, and computed at the position of the vertex (x^ = 0), is zero: 

Y, a nf*toM + (w T w)^n lcad ^) = o . (22) 

(3 

We have used the obvious notation: f^(^) = ^^(x^ = 0). The second term is the con- 
tribution of a lead, if one is plugged at vertex \i [due to the definition of W, we recall that 
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(W T W)^ = 1 if a lead issues from \i and otherwise]. The second useful property is obtained 
by integration of ( |2l| ) on the bond 

dx \^ W) (x)\ 2 = -SltfQi) - n^03) . (23) 

We now consider the stationary scattering state ip( a '(x) of energy E = k 2 , associated to 
the injection of a plane wave from the lead connected at vertex a. The construction of these 
eigenstates is briefly recalled in appendix ^ (see From the expression ( |47| ) of the wave 

function on the lead we see that 

n&>) = -^U^w ~ k ^ a + ^ { ~ 6 ^ + ■ (24) 

We now compute the integral of \if}^ a \x)\ 2 on the graph (the "graph" refers to internal bonds): 



L 



dx\^ a \x)\ 2 = £ ! dx\i>^ p) {x)\ 2 =-J2 , (25) 



Gra P h J arc jiv 



where we have used (|23|). The summation Y2(u/3) * s over the ^ bonds of the graph whereas 
the last summation runs over the IB internal arcs. We see that the contributions from the 
arcs issuing from an internal vertex vanish due to (|22|). The contributions of the internal arcs 
issuing from connected vertices can be replaced by the contributions of external leads due to 
(l22l). Therefore we obtain: 

f dx\^(x)\ 2 = ^0{^>) = ~2ik (st^|) - ^(S Qa - E*J (26) 

J Graph „ \ / ctot 

where the sum over /i is obviously over the L connected vertices. 

In order to associate a measure dE to the stationary scattering states, we change the normal- 
ization. The scattering states (47,50) are related to the new ones by: ?/^(x) = ^-j- ip^(x). If 
we sum the contributions (p6|) of the L stationary scattering states of energy E, we obtain: 




which generalizes the Smith relation ( |i~8|) to the case of graphs. The term — iTr {S^^j is the 
time delay. 

To compute the Friedel phase of a graph appearing in the above relation, it is useful to note 
thatfl 

det(W?W-M) 
detS - (_1) det(W?W + M)- (28) 



3 The proof is easily achieved by considering the graph Q' related to the original graph Q by attaching to each of 
the V — L internal vertices of Q (labelled for convenience with prime indices: a ,. . . ) a lead with tunable coupling. 
If these couplings are swichted off (w a i — > 0), the V x V scattering matrix E' of Q' is block diagonal with a L x L 
block being the scattering matrix E of Q, the other (V — L) x (V — L) block corresponding to the additional leads 
being —1. Let us now compute detE': for finite couplings w a i , the matrix W' describing the coupling of Q' to 
the V leads is square and possesses an inverse. It follows from (|l|) that E' = W' _1 (W' 2 - M)(W' 2 + M)~ l W'. 
Then detE' = ^"t(w ,: ^+Aj) ' ^ ^ e cou pli n g s to the additional leads now vanish, w a i — > 0, we have detE' = 
(-l) v_i detE and W' 2 = W T W. Qed. 
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For one channel (one lead) we have S = e 1<s , therefore we get 

L^°w 2 =H^ + ^ s )' <29) 

which would be the relation (|l8|) obtained by Smith if the graph would reduce to a line (one- 
dimensional case with one channel). The different sign of the second term is only a matter of 
definition of the phase shift 5, in the one-dimensional case, and 5q, in the I = channel of the 
three-dimensional case, which differ by it. 

Case w a 7^ 1 

When we introduce arbitrary couplings between the leads and the graph, the application of 
formula (^7|) means that we are also taking into account integral over the bonds on which are 
the barriers characterized by w^s (see Jl) where the introduction of these parameters is explained 
in detail). 

4 Violation of the Friedel sum rule for certain graphs 

The idea of the FSR is to count the states in the scattering region by studying the Friedel phase. 
We have seen that the Smith formula (|ilj|) relates the LDoS integrated over the scattering 
region to the Friedel phase and we have found its generalization ( p7| ) for graphs. We call 
M(E) = J^qq dE' J Graph dx p(x; E') the integrated density of states (IDoS) of the graph. If we 
are not interested in the details of the spectrum but only in the Weyl term of the IDoS of the 
scattering region, and if we believe the FSR, the relation ( |l9| ) shows that M^ ey \(E) ~ ^8 ■'(E) 
up to some oscillatory part. As a matter of fact this is not always true for graphs and we will 
now give several examples where A/weylC^O is n °t given by the dominant contribution of 5^(E). 

All the examples we are going to consider are free graphs, with V(x) = 0, but the ideas that 
will come out are not specific to free graphs. We recall that in the absence of a potential we 
have^: Afw ey i(E) = where C = | Yla 8 a a(3^af3 is the total length of the graph. As we will see, 
one of the reasons of the violation of the FSR is the occurence of degeneracies in the spectrum. 
Graphs with symmetries present a lot of degeneracies. This is why it is interesting to start by 
studying the complete graph Ky, which is the most symmetric simply connected graph with V 
vertices. 

4.1 The complete graph Kv connected to one lead. 

The graph Ky is made of V vertices, each being connected to the other ones by bonds of same 
length I. The matrix M takes a simple form (see equation ( |64| ) in appendix ^). 



Figure 1: Complete graph K§ connected to one lead. 
4 The Weyl term appears in the trace formula originally derived by Roth 0] (see also [0, W). 
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The graph is connected to one lead (figure ||) and the scattering matrix is characterized by 
a unique phase: £ = e 1 ^ . Using (|l3|) we recover after a little bit of algebra the expression 
formula (119)]: 

COS ki + COS 09 — 1 . , „ , . . . 

cotg^/2) = cos ip -— — t cotg (k£ 2) , 30 

cos ki + cos ip 

where k = \^E and cos ip = p^ry- This expression shows that 8$ (k 2 ) = 3k£ + (fluct.), where 
(fluct.) represents a fluctuating term of order n, whereas the Weyl part of the IDoS is J\fw ey i(E) = 
B^k _ V(V-l)k£ ^ w j 1 j c ] 1 c i ear jy shows the discrepancy between M{E) and • 

A more detailed analysis of the position of the resonance peaks of shows that the Friedel 
phase does not measure the degeneracies of the energies of the isolated graph (see appendix 
|c] where the spectrum of Ky is recalled) and moreover even misses some energies: there is no 
resonance peak at fc2+4n- 

To understand in more general terms the origin of the failure of the FSR when only one lead 
is plugged on the graph, we consider the simple case of a graph with no potential (V(x) = 
and X a = 0) connected to only one external lead. The formula *M 

relates the phase shift 8$ (E) to the ratio of two spectral determinants. On the one hand 
Star. (7) is the spectral determinant det(— D^+7) calculated with a Dirichlet boundary condition 
(A ao = 00) at the vertex ao where the lead is plugged in, and Neumann boundary conditions at 
all other vertices (X a = 0). On the other hand 5N eu . (7) is calculated with Neumann boundary 
conditions at all vertices. The sum rule means that each state in the isolated graph is associated 
to a jump of 2tt of the phase 6*. Due to (31) we see that a jump of 2tt occurs when the expression 



\3l\ ) diverges. Then we identify two reasons why the FSR fails: (i) if the spectrum of the graph 



is degenerate and (ii) if Sbir.(7) vanishes for the same energy as <SNeu.(7)> then (31) diverges a 
number of times which is not related to the number of states in the graph. 



4.2 The complete graph K v connected to V leads. 

To convince ourself that the breakdown of the Krein-Friedel relation is not specific to graphs 
connected to one lead only, we consider now the case where Ky is attached to V leads connected 
to each vertices (figure |2|) . 




Figure 2: Complete graph Kj connected to 7 leads. The small boxes represent the couplings 
characterized by the parameter w fij. 
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If all the vertices of the graph are connected to leads (L = V), the matrix W is square. The 
determinant ( |28| ) is: e lSf ^ = dct(w^+M) • ^ or simplicity we consider the case of equal couplings: 
w\ = W2 = ■ ■ ■ = wy = w. Using (|6i]) we get: 

(fc 2 ) / wv sin(M/2) — \w 2 cos (pcos(k£/2) / cos k£ + cos ip + ku 2 cos ip sin H 
sin(M/2) + iu; 2 cos 99 cos(k£/2) \cos kl + cos 99 — iw 2 cos 99 sin kl 

This expression shows that S f {k 2 ) = (2V - l)k£ + (fluct.) which disa grees once again with 

It is interesting to provide a more detailed analysis by studying the behaviour of the Friedel 
phase in the neighbourhood of the energies of the graph (the spectrum is recalled in appendix 
|c|). We consider the limit w — ► for which the resonance profile of t™- emerges clearly. 

• Near the first energy level (for k ~ ki) we see from (|32|) that e 1<5/ tx f fE^tpjaf^ 1 with 
Ak±£ = w 2 /(V — 1) . The exponent is the degeneracy of the level which means that the resonance 
peak of ^ has the correct spectral weight and counts correctly the V — 1 states. 

• In the neighbourhood of the second energy level (for k ~ k 2 ) % is flat: 5 f is not sensitive to 
the presence of states at this energy. 

• The situation at k ~ k% is the same as the one at k ~ k\ (with Ak^ = Ak\) 

• At k ~ k^ we have e 1<5/ cx with Ak^ = 2Ak\: the Friedel phase misses the degeneracy. 

One may now ask the question why the Friedel phase misses some states sometimes and 
sometimes not ? To answer this question we can study the structure of the wave functions of the 
isolated graph (see appendix 0). Whereas the wave function is finite at the nodes at energies 
k\+2m where the Friedel phase is sensitive to the degeneracy, all the V(V — 3)/2 degenerate wave 
functions vanish at all the nodes at energies &;2+4m as well as at energies /c4 + 4 m , which means 
that the wave sent from the lead does not enter the graph at those energies. 




4.3 The ring connected with two leads. 

To understand better the remark that closed the previous subsection we consider next a simpler 
case: a ring connected to two leads (figure |3|). The arms of the ring are of length l a and lb with 
I = l a + lb- 

Let us first recall the spectrum of the isolated ring of perimeter I threatened by a flux 9: 
the energies are E m (9) = (^) (m — ^) associated to wave functions <p m (x) = ^Q 2mmx / 1 f or 
m G Z. If = the states ip m and (p- m are degenerate. Therefore we can introduce in this case 
a different basis: a symmetric function <p£ = y "^|~ n = \/2/l cos(2irnx/l) with n £ N, and an 

antisymmetric one: (p~ = yn 7^~" = y/2/l sm(2imx/l) with n € N*. 




h 



Figure 3: The ring connected to two leads and threatened by a flux 6. The two arms have 
lengths l a and The parameters w\ and W2 allow to tune the coupling of the ring. 
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We now consider the graph when it is coupled to two leads with coupling parameters w\ and W2 

det(W 2 -M) 
&et(W 2 +M) 

iS f 2(cos 9 — cos kl) + w\ sin kl a sin klb — i(w{ + w^) sin kl 
2(cos 9 — cos kl) + w\ w\ sin kl a sin Mb + i(u> 2 + u^) sin kl 



(see figure H). The Friedel phase e 1<5/ = det(iy 2 +M) is: 

iS f 2(cos 6* — cos kl) + u?^ t/?2 sin kl a 0111 rui.£, — i v u -'i ~r~ uy 2y' di±a >f 

e = n zzzri^~r~2~^TznrTr~. rm , :^..2 , „..2\ m n (33) 



(det(PF 2 + M) was calculated in If # 7^ the spectrum is non degenerate and <F counts 
correctly the states. Now we focus on the degenerate case 9 = for which we have: 

tan(tf//2) = (w} + v%)*nkl 

4 sin (kl/2) + w(w2 sin kl a sin kit, 

Each interval of width A/c = 27r// contains 2 states of the ring. Let us now examine under what 
condition ( |34[ ) counts correctly these states. We can identify the position of the resonances with 
the value of k for which the denominator of ( |34|) vanishes^]. Then we distinguish two different 
cases: 

• l a /l is an irrational number. Then the denominator of (|34]) vanishes twice per interval k G 
[2mir/l; 2{m + l)ir/l[, which means that 6? counts the correct number of states. 

• l a /l is a rational number: l a /l = with (p, q) G N 2 . If (q — p)(m + 1) is an integer multiple 
of q, the denominator vanishes only once in [2rrnr/l;2(m + 1)tt/Z[. The intervals for which 
(g — p)(m + 1) is an integer multiple of q are those in which one of the two degenerate wave 
functions and vanishes on the vertices where the lead are plugged in. 

4.4 The ring connected with one lead. Why the l.h.s of (|27]) can not always 
be identify with the DoS ? 

We have given a general argument to explain how the degeneracies of the spectrum lead to a 



failure of the FSR, however it is surprising that the quantity in the left hand side of (27) can 
not always be identified with the DoS of the graph since the sum runs over the complete set 
of L stationary scattering states of energy E. This point needs a clarification that we will give 
now by studying again the case of the ring. 




Figure 4: 

We consider the ring of figure |3] coupled to a lead and construct the stationary scattering 
states (47j50|). The wave function on the lead is 

^ ad (x) = ( e~ ikx + e ikx+iS ) , (36) 
\/4irk \ ' 

The denominator is of the form: 

f a (x) = sin 2 (x/2) + bsin(ax) sin((l - o)a;) = (1 + 6) sin 2 (a;/2) - 6sin 2 ((l/2 - a)x) , (35) 

with a g]0; l/2[ and b g]0; +00 [. We are interested in the number of zeros of f a (x) in the interval [2mn; 2(m + l)n[. 
We distinguish two cases: 

• a is not a rational number (a £ Q). Since / a (2m7r) < Vm g N and the amplitude of the first positive term in 
the r.h.s of ( ^Bj ) is larger than the second, it follows that f a (x) — has exactly two solutions in [2rmv; 2(m + l)7r[. 

• a G Q : we write a = 2- where (p, q) £ N 2 with p < q. We have / a (2m7r) = — &sin 2 ((g —p)mir/q) ^ 0. 

If (g — p)(m + 1) is not an integer multiple of q the interval [2m7r; 2(m + 1)tt[ contains two solutions of f a (x) = 0. 
If (q — p)(m + 1) = rq with r £ N the interval contains only one solution of f a (x) — 0. 
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where the phase shift 5 is H 



, „ , . w 2 sin kl , . 

= 2( c osfl-cos«) ■ < 37 > 



On the ring (arc a) the wave function reads: 

V'Cal (x) = ibx — — — ( si n k(l — x) + e _1 sinkx) (38) 
v ; sin/a V / 

where 

ft-' ,. 51=^(0) (39) 

V vrA; if sin fct + 2i(cos fct — cos v) w 

is the wave function at the node. 

• 0^0. We study the limit w — > of small coupling for which we expect to recover some 
features of the isolated ring. In this case |^i| presents sharp peaks at the positions of the energies 
of the isolated ring (given by k±l = ±9 + 2mr). These resonance contributions will eventually 
give the main contributions to (|27|). Let us express the wave function in the ring for k in the 



neighbourhood of a resonance. Expressions (38,33) give: 



fc~fe± V-xk k-k n + iw z /2l 

Up to a normalization, we recover the wave functions of the isolated ring recalled in the previous 
subsection. By integration in the ring we obtain: 

f dx^ {a) (x)\ 2 g i- - ^ — » i-f (fc - kt) = 5{E - \kt) 2 ) , (41) 

Jo fc~fc± 2k n (k — kn) + (w 2 /2l) z w-*o 2k 

which is the correct DoS of the isolated ring. 

• = 0. If we now consider the degenerate case of zero flux, the resonance peaks are in 
k n = 2nir/l. The wave function in the ring near the resonance is 

W)(*) * -7= , r'l- 2/7 cos(2nW0- ( 42 ) 
fc~fc„ V7rA; k-k n + vw l /l 

In this case the scattering state only reproduces the symmetric wave function (fn(x) of the 
isolated ring. It is now clear that the integration can not give the DoS of the isolated ring: 
indeed, 



/ dx\^ (a) (x)\ 2 — —r - 



u; 2 // 1 



—S(k - k n ) (43) 



(k - k n ) 2 + (w 2 /l) 2 w -, 2k 
misses the degeneracy 2 of the eigenstates. 

In the cases studied above the stationary scattering states, computed in the limit of a graph 
weakly coupled to the leads (w a — * 0), do not reproduce all the wave functions of the isolated 
graph and the l.h.s. of the formula can not be identified with the DoS of the graph. In this 
sense the scattering states do not form a complete basis to describe the Hilbert space of the 
graph. 

5 Remark on free graphs connected to one lead 

Note that for a free graph, i.e. with no potential and no magnetic flux, perfectly connected 
(w = 1) to only one lead, we can prove with pi]) that 



5 f (k 2 ) ~ 2kC (44) 
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where C = 5 ^ Q a a a pl a p is the total length of the graph. The proof is obtained by analyzing 
the behaviours of the two spectral determinants at small energy R. If the FSR fails, the total 
length C of the free graph is not encoded in the Weyl part of 5^(k 2 ) but however £ appears in the 
low energy behaviour of the phase shift. For example if we consider the complete graph (figure 
[l]), we have shown that S{^ eyl = "Ski and we can check on ( |30|) that the low energy behaviour 
(k -» 0) is 5 f ~ V(V - l)k£ = 2kC. 



6 Discussion 

We have shown that the well-known Friedel sum rule (FSR), a state counting method from the 
scattering properties, may be violated for certain graphs having degenerate spectrum. This has 
been demonstrated already at the level of the Weyl term of the DoS: we have studied several 
examples where a discrepancy occurs between the Weyl term of the Dos of the isolated graph 
and the Weyl part of the derivative of the Friedel phase ^- 4r|r = ^ Tr { £t } . 

A way to understand the origin of the failure of the FSR is to compare the quantities involved 
in the DoS and in the Friedel phase. The density of states of an isolated graph can be obtained 
from the spectral determinant S(j) = IJn(7 + E n ) by doing the substitution 7 — > —E — i0 + . 
The spectral determinant is proportional to the determinant of the matrix M introduced above: 
5(7) oc det M(j) |], |34], [7], 10, |T^| . Adding a small imaginary part to the spectral parameter 7 
produces a resonance structure in <9 7 In S(-y), each peak having a weight equal to the degeneracy 
of the state. If we now consider the Friedel phase we note that the width of the resonances are 
obtained by adding to the anti-hermitian matrix M(—E) an hermitian matrix W T W: the Friedel 
phase involves det(M(—E) + W T W). Comparing this latter determinant with the determinant 
det M{—E — i0 + ) involved in the density of states, it is not surprising that the Friedel phase 
does not produce the correct spectral weights since the ways the energies (zero of determinant) 
acquire an imaginary part is different in the two cases. 

Another way to understand the failure of the FSR for graphs is the following. For a problem 
invariant under rotations in a d-dimensional space, the essence of the FSR is to count the number 
of nodes of the wave function in the angular channel crossing a {d — l)-dimensional sphere at 
infinity when the energy is varied. The number of states coincides with the number of nodes, 
that is with the number of jumps of n of the phase shift rji(E) of the partial wave of orbital 
momentum I. On the other hand, a graph is connected to the external only through leads plugged 
at vertices. In a sense the Friedel phase counts the number of nodes of the wave function ipk(x) 
that reach those vertices by varying k. The failure of the FSR is caused when several nodes of 
ipk{x) reach at the same energy the same vertex from different bonds issuing from this vertex 
(we can easily convinced ourselves of this remark by considering the ring connected to one lead 
studied above). 

We can also provide a clear picture of the problem within the arc formulation introduced in 
|Q and recalled in section |l]. In the arc formulation, the wave function is described by a set of 
amplitudes. Each arc i is associated with a couple Ai,B{. We gather the internal amplitudes 
in a vector A mt , the external in a vector A ext and all amplitudes in a vector A. The internal 
amplitudes of the graph are related through the bond scattering matrix : A mt = RB mt . All 
amplitudes are also related with each other by the vertex scattering matrix : B = Q A. If we 
eliminate B mt we get : 

Q T A ext = (i?t - Q int ) A int (45) 
B cxt = Q A int + Q cxt A cxt . (46) 

In general det(R) — Q mt ) 7^ whatever k is and at all energies of the continuous spectrum, 
the stationary scattering states are the only solutions of the Schrodinger equation on the graph. 
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However, for certain graphs (in particular for those examined above), there exists a discrete 
set of energies in the continuous spectrum for which det(R) — Q mt ) = 0. This means that 
additionaly to the scattering states, we can construct at those particular energies solutions 
such that ^4 ext = B cxt = while the internal amplitudes satisfy (1$ — Q mt )^4 int = and 
QA mt = 0. These two last equations describe a solution localized in the graph and that does 
not communicate with the leads. The stationary scattering states give the solutions of the 
Schrodinger equation for the continuous spectrum apart for a discrete set of energies where 
some additional states are localized in the graph and thus are not probed by scattering leading 
to the failure of the state counting method from the scattering. 

The study of the various examples of section ||| leads us to make the following conjecture 
for the ability of 6* to count the states (at least at the level of the Weyl term): if there are 
degenerate energies of degeneracies d n , the Friedel phase 5* counts correctly the states of the 
system if L ^ d n leads are plugged at vertices in such a way that the wave function can not 
vanish at the positions of all these vertices at the same time. 
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A The stationary scattering states 

In this appendix we recall briefly how the stationary scattering states are constructed 0]. We 
consider the stationary scattering state ip^(x) of energy E = k 2 which describes a plane wave 
entering the graph from the lead connected at vertex a and being scattered by the graph into 
all leads. On the lead connected to vertex fj,, the wave function is: 

Wl>) = ^e-^ + £^e ifc *, (47) 

where x E [0; +oo[. The wave function on the internal bond ((10) of the graph, is related to the 
two linearly independent solutions f lx p(x^ L p) and fp^x^p) of the differential equation 

(-d^ + Vw)M - k 2 ) fM = (48) 

for x E [0; l^p]. The two solutions f^p and fp^ satisfy the following boundary conditions at the 
edges of the interval: 

(//*(m) = i and = o 

\ = o \ um = 1 

where /(//) = f{x jJb p = 0) and f{(5) = f{x lx p = l^p). The stationary scattering state on the bond 
((10) is: 

= ^ ] fi*M + < } fp*M ( 5 °) 

which already satisfies the continuity condition (§). The relation between the functions f^p 
and fp^ and the reflection and transmission coefficient characterizing the potential on the bond 
is established by computing the derivatives of f^p and /«„ at the boundaries of the interval. 
Imposing the "current conservation" (||) then leads to the expressions (13,12,|ll]) that permit a 
systematic construction of the scattering matrix O. 
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B Friedel sum rule and Smith relation in one dimension 



The one-dimensional case can be considered as a graph with one bond and two vertices and can 
therefore be described with the formalism presented in O and this paper. The FSR (|l~6|) has been 



demonstrated in general terms in [15], however it is interesting to give a rapid demonstration 



that follows the lines of the original one in three dimensions [14, 16, 20, ^] ; note also that it 
has been demonstrated in that the Friedel phase in a one-dimensional situation is related 
to the phase of the transmission amplitude (see also |30|). We consider the one-dimensional 
hamiltonian — d 2 + V(x) with i£l with a potential V{x) being concentrated in some region of 
the space. We start by describing several possible basis of eigenstates caracterizing the scattering 
problem. 

• The stationary scattering states of energy E = k 2 related to the scattering matrix 

H't'?) (5l) 

are the state <p( L \x) associated with a plane wave coming from the left and ip( R '(x) for an 
incoming wave from the right. The asymtotic behaviours of the left stationary scattering state 
are ^ L \x) = e lkx + re~ lkx for x — > — oo and tp( L >(x) = te lkx for x — > +oo. The state ^ R \x) 
involves similarly the coefficients f' and t' . Note that those states can be introduced even if the 
potential is not concentrated in a finite interval provided that it decreases sufficiently rapidly at 
infinity. 

• If the potential has a support [xi;x2] we introduce the stationary scattering states if)^ L '(x) 
and ip( R \x) related to the scattering matrix S: the left stationary scattering state behaves like 
ip( L \x) = Q lk ( x ~ x i) _|_ re -ik{x-x 1 ) £ or x ^ Xl anc } ^( L )(x) = te lk( - x ~ X2 ^ for x ^ X2- A similar 
expression for ip( R \x) involves the coefficient r' and t'\ i[)( R \x) = tfe- lk ( x - x i) f or x ^ x\ and 
^ R \x) = e~ lk ( x ~ X2 ^ + r ' Q lk ( x ~ x 2) f or x ^ X2- The reflexions and transmissions defined in this 
way are naturally involved in transfer matrices, which makes one of the interest of this definition. 

Comparing the two sets of eigenstates it is clear that tp <yL \x) = e lkxi ^ L ^ (x) and ^ R \x) = 
e~ lkX2 ijj( R \x). The relations between the coefficients of the two scattering matrices S and S 
are then: r = f e~ 2lfea;i , r' = f'e 2lfcx2 , t = te lk( - X2 ~ Xl ^ and t' = p e lk ( X2 ~ Xl \ The relation between 
matrices reads: £ = UYM with U = diag(e~ lA:Xl , e lkX2 ). 

• To derive the FSR we introduce the two eigenstates of energy E = k 2 labelled by the index 
a = 1,2: 

^ a {x) = [<V+ 9(x) + a a - 0{-x)} sm(k\x\ + f] a (k 2 ) + ir/2) for \x\ ->■ oo , (52) 

where 6{x) is the Heaviside function. If the potential is symmetric the two amplitudes a a<+ 
and a a - are equal in modulus and a = 1,2 labels the symmetric and antisymmetric states. 
Let us establish the relation with the 2x2 scattering matrix S. We look for the relation 
between this basis of eigenstates and the first basis introduced above: let us write ^ a (x) = 
(p( R \x) + C <p( L \x). Comparing their behaviours at x — > +oo we get f' + Ct = e 2l??CT and at 
x — > — oo: t' + Cr = C e 2l??CT . These two equations show that e 4ir?CT — (f + f')e 2l?7CT — tt' + rr' = 0. 
In other terms: 

det(S-e 2i ^) = 0. (53) 

Therefore e 2ir?1 and e 2ir?2 are the two eigenvalues of the scattering matrix S. To finish the proof 
of the FSR we consider that the system is in a large interval [— R; +R] and impose Dirichlet 
boundary conditions. The quantification condition for ^ a (x) reads k n R + r\ a {k^) + 7r/2 = nir. 
We introduce 5k n = k n+ \ — k n , therefore in the limit R — > oo: 

1_^R + 1 d Va (k 2 n ) 

Sk n 7T 7T dk n 
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which is the density of modes in the channel a. The term ~ is the density of modes in the absence 
of the potential: . In the limit R — > oo the difference of densities of modes ^ — remains 

finite. It follows that 

f + ™ dx [p(x; E) - p (x; E)} = I £ (55) 



<T=1,2 



where p(x; E) = (x \5(E — H) \ x) is the LDoS and po(x; E) the LDoS in the absence of the poten- 
tial. Due to (|53|) this equation can be rewritten f^^dx [p(x;E) — po(x;E)] = ^ Tr j ^ | . 
Qed. 

Next we would like to apply both the FSR and the Smith formula on a simple example. 
We now consider a potential with support [0; C] that vanishes elsewhere, a situation where it is 
meaningful to introduce £ (instead of £). 
(i) The Smith formula (|27| ) gives the DoS of the interval [0; C] 



J dxp(x;E) = dx(|^ L) (x)| 2 + |^ ) (x)| 2 ) = i-(Tr|st^} + i LTr{s-St}) . 

(56) 

The two terms correspond to left (L) and right (R) stationary scattering states, which form a 
complete basis of eigenstates in one dimension. Note that this relation has also been given in 



[28] for the one-dimensional case. 



(ii) On the other hand the FSR 

measures the variation of the DoS of the infinite line due to the presence of the potential in the 
interval [0;£]. In particular ( |57| ) is sensitive to the effect of the potential on the wave function 
at infinity whereas (^) is a local quantity. 

As an illustration, let us consider the extremely simple case of a potential X5(x) on a line. 
The corresponding graph is a vertex (V = 1, B = and L = 2). Formulae (^J^,||) give the 
scattering matrix 

2 fl ^ 1. (58) 



2 + \\/k V 1 1 
It is easy to check that (^) therefore vanishes 



(59) 



which is not surprising since the "graph" is only a point (the support of the potential is an 
interval of measure 0). On the other hand the variation of the DoS of the infinite line can be 
computed either with the exact Green's function, known for this potential, or with (|57|). The 
Green's function gives: 



/+oo 1 A 1 

dx [p(x; E) - P0 (x; E)] = 8(-\)5(E + A 2 /4) - -6(E) + d ( E ) E + X 2 /A » 



(60) 



where 9(E) is the Heaviside function. The first term is the contribution of the bound state (that 
exists if A < 0). We can check that the total number of states is conserved: 
f-m d-E 1 f^°° dx [p(x; E) — po(x; E)] = 0. The Friedel phase, obtained from the above scattering 
matrix is 8* (E) = — ilndet S = 2 arctan(2/c/A) and we can therefore recover the expression (|60|) 
using the FSR @). 
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C The spectrum of the complete graph K v 



We give here the spectrum of the complete graph, which is made of V vertices all connected 
among each other with bonds of same lengths £. The spectrum is easily extracted from the 
spectral determinant S(—k 2 ) = []™ =0 (£ m — k 2 ) which has been computed in |7|: 



S{-k 2 ) oc 



sin k£ 
k 



V(V-3) 



sm 2 {k£/2) (cos k£ + cos ip) 



v-i 



(61) 



up to some inessential numerical factor. The parameters X a that characterize the boundary 
condition f^) are put to zero here. We have introduced cos (p = yry- The energies E m and the 
corresponding deg ener&cies dm £ire given in the following table: 





dm 


k = 


1 


h = ^ 


V-l 


k 2 = \ 


V(V-3) 
2 


ft 3 — £ 


V-l 


k„ - 22L 
K 4 - e 


2 + v ^ 



It is obvious from the expression of the spectral determinant that the spectrum is periodic in k 
of period 2ir/£, that is: /c m +4 = k m + 2ir/£ and d m +i = d m (for m > 0). 

We next consider the corresponding eigenfunctions. The eigenfunction on the bond (otf3) is 
given by (M): 



VV/3)( X ) 



1 



(ijj a sin k(l a/ 3 — x) + ip/3 sin kx) 



sin kl a g 

Imposing the conditions (||) leads to the V equations: 

P 

For a free graph M is given by (|l4|) P|. For the complete graph we have: 

i 

{S a j3[V — 1) cos k£ — a n 



(62) 



(63) 



M, 



a/3{ 



-k Z 



sin k£ 



la/3) 



where the adjacency matrix is a a p = 1 — 5 a p. 

• Zero Mode. The wave function is constant on the graph (x) 
is the total length of the graph. 

• k 



1 / VZ. where C 



(64) 



V(V-1) . 



The equation (|63|) has V — 1 

S a ,i ~ <Wfi, up to a 



fci. All the matrix elements of M are equal: M a p = — 

(1 i) 

solutions labelled by j = 1, 2, • • • , V — 1. A possible basis is: ipa 
normalization (this basis is not orthogonal). 

• fc = fc 2 - At this energy the matrix M is divergent and equation (^) can not give the 
eigenstates. They are obtained by considering the equation (1 — RQ)A = where R and 
Q are the matrices given by (^j6j,RM. A is the vector gathering the 2B amplitudes of the 
wave function (one for each arc) The system (1 — RQ)A = has V(V — 3)/2 solu- 

tions at k = k.2- To have an idea of the structure of the solution let us consider K4. We 
label by 1,2,3 and 4 the nodes, and we bring together the 6 components on the 6 bonds in 
a vector V(x) = (V>(12)0»0, ^(13) 0*0^(14) 0*0 > ^(23) 0*0 > ^(24)0*0^(34) 0*0)- We have for the first 
state, labelled (2,1): ^I^ 2 ' 1 ) = (1,0,-1,-1,0,1) x sin(irx/£), and for the second eigenstate: 
^(2,2) = ^ 1; _ x _ x 1; j x sin ( nx /£y 
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ip(2,l) ip(2,2) 




Figure 5: The two eigenstates of K4 of energy k = &2- The dotted lines are the bonds and the 
large dots the vertices (labelled 1,2,3 and 4). The large continuous lines are put where the wave 
function is positive and the large dashed lines where it is negative. On the other bonds the wave 
function vanishes. 

• fc = fc 3 . The matrix M is the opposite as the one computed at k = k\ and the wave functions 
on the nodes have the same value as for this latter energy. 

• k = &4. The same problem occurs as for k = /c2- The system (1 — RQ)A = has V(V — 
3)/2 + 2 solutions corresponding to wave fonctions vanishing at all the nodes. Again we consider 
the graph K4 and give the four degenerate states: ^I^ 4 ' 1 ) = (0,0,0,-1,1,-1) x sin(2-7rx/^), 
*(4,2) = (0,-1,1,0,0,-1) x sin(27rx/£), f( 4 > 3 ) = (-1,0,1,0,-1,0) x sin(27ra;/^) and ^( 4 < 4 ) = 
(-1,1,0,-1,0,0) x sin(27rx/£). 



\p(4,l) \p(4,2) \p(4,3) \p(4,4) 




Figure 6: The four eigenstates of energy k = k^. 

• k = k m> 4. The spectrum is periodic in k with period 2ir/£. Then the value of the wave 
functions at the nodes is the same at k m and k m+ 4, only the number of oscillations on the bonds 
change. 
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